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Abstract: We study gravity in backgrounds that are smooth generahzations of the Randall- 
Sundrum model, with and without scalar fields. These generalizations include three-branes in 
higher dimensional spaces which are not necessarily Anti-de Sitter far from the branes, intersecting 
brane configurations and configurations involving negative tension branes. We show that under 
certain mild assumptions there is a universal equation for the gravitational fiuctuations. We study 
both the graviton ground state and the continuum of Kaluza-Klein modes and we find that the 
four- dimensional gravitational mode is localized precisely when the effects of the continuum modes 
decouple at distances larger than the fundamental Planck scale. The decoupling is contingent only 
on the long-range behaviour of the metric from the brane and we find a universal form for the 
corrections to Newton's Law. We also comment on the possible contribution of resonant modes. 
Given this, we find general classes of metrics which maintain localized four- dimensional gravity. 
We find that three-brane metrics in five dimensions can arise from a single scalar field source, and 
we rederive the BPS type conditions without any a priori assumptions regarding the form of the 
scalar potential. We also show that a single scalar field cannot produce conformally-flat locally 
intersecting brane configurations or a p-brane in greater than {p + 2)-dimensions. 



*J. Robert Oppenheimer Fellow 



1. Introduction 



The proposal of Randall and Sundrum (RS) [1,2] to localize gravity in the vicinity of a brane with 
non-vanishing tension in anti-de Sitter (AdS) space has recently attracted enormous attention 
(see, for example [3,4], for previous relevant work, [5-24], for more recent generalizations, [25-30], 
for work on smooth brane scenarios, [31-37], for embeddings in string theory and supergravity, 
[38-47], for the general relativity aspects and finally [48-52], for cosmological and phenomcnological 
aspects). RS found that in a setup with a single brane, a negative bulk cosmological constant and 
a single large extra dimension (with the cosmological constant and brane tension tuned such that 
the effective four-dimensional cosmological constant vanishes) the solution to Einstein's equation 
results in a single graviton zero mode, which is a consequence of the unbroken four-dimensional 
Poincare invariance, and a continuum of Kaluza-Klein (KK) modes. Normally the presence of 
these continuum modes would render a setup like this unrealistic due to the large deviation from 
Newton's Law the low energy continuum modes tend to induce. However, RS found that due to 
the suppression of the wavefunctions of the continuum modes close to the brane, their contribution 
to the Newton potential is highly suppressed, and therefore a realistic model with uncompactified 
extra dimensions could be built. This model has been generalized in [5] to models with intersecting 
branes with more than one uncompactified extra dimension, and also to include brane junctions 
[6,7]. 

The branes in the RS setup and its generalizations mentioned above are included as static 
point-like external sources in the extra dimensions, with no dynamics to produce them. As was 
done in [26-29,33], one can find solutions to Einstein's equation coupled to a single scalar field, 
where the scalar creates a domain wall — a "thick brane" — while the metric away from the brane 
asymptotes to a slice of AdSs. Such domain wall solutions are obtained if the scalar potential 
originates from a superpotential (although as recently discussed in [37] this does not necessarily 
imply that the theory is embeddable into a five-dimensional supergravity theory) . In this case the 
solutions found in [26-28] originate from a BPS equation. These domain walls were first found 
in [25]. It has been shown in [26-28] that, just like for the case of the infinitely thin branes of RS, 
there is a single normalizable graviton bound-state with zero energy. A particularly nice example 
of this sort has been recently worked out in detail in [30]. Similar BPS equations for intersecting 
domain walls in more than one extra dimension were found in [8, 9] ; however, no explicit solutions 
to these equations are known yet. 

In this paper we study generic properties of localized gravity on thick branes. In the first 
part of the paper we consider thick branes in one extra dimension and then generalize to an 
arbitrary number of extra dimensions. Instead of starting with a coupled gravity-scalar system, as 
in [27,28], we "smear" the RS solution and its generalizations in such a way that the non-dynamical 
source terms correspond to a smeared (thick) brane in the background of a slowly varying negative 
bulk cosmological constant. We examine the spectrum of graviton modes and find necessary 
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and sufficient conditions for such backgrounds to locafize gravity on the branes. Besides general 
arguments about the ground state (some of which appear in [27,30]), we also examine the behavior 
of the continuum modes. For a generic study, we use the WKB approximation for the "volcano- 
type potential", which hints that when the metric falls off slowly enough from the brane the soft 
KK modes are suppressed sufficiently inside the brane so that their corrections to the Newton's 
Law are negligible. In a more restricted set of generalizations of the RS solution we calculate this 
suppression more rigorously, and find qualitative agreement with the WKB result. This leads to 
necessary conditions of the quantum mechanical potential, and hence the background metric, in 
order for the KK modes to decouple. We find that the potential at large distances must fall off no 
faster than in the asymptotically AdS case in order for the KK modes to make a small contribution 
to Newton's Law. This requirement is equivalent to demanding normalizability of the ground state 
graviton wavefunction. We also comment on the possible contributions of "quasi-bound-states" — 
resonant modes in the continuum spectrum whose wavefunctions are not suppressed at the location 
of the brane — and study their significance in a toy model. We next show how to generalize our 
results to situations in more than five dimensions. These scenarios could describe, for example, 
three-branes in more than five dimensions or higher-dimensional intersecting branes with a four- 
dimensional intersection. In the latter case, the thick brane background could be given by an 
appropriate smearing of the intersecting brane solution of [5] , and again we find conditions on the 
long-distance behavior of the background metric in order for there to be localized gravity on the 
brane intersection. 

We also study the relevance of background fields that create the branes. We find that the stress- 
tensor source terms for a general smearing of the RS solution can be obtained from a single scalar 
field, and we rederive the same BPS- type equations for this scalar field as [26-28]. This provides 
a particularly simple derivation of the BPS equations without an a priori assumption about the 
form of the scalar potential, and also emphasizes that this is the most general solution with a single 
scalar field. In the case of branes in higher dimensions the situation is more complicated. Contrary 
to the case of one extra dimension, we find that it is impossible to generate the desired background 
metric or sources of the stress tensor with a single scalar field. Nevertheless, the properties of the 
graviton in such backgrounds are studied in the same way as for the case of one extra dimension. 



2. Backgrounds with four- dimensional Poincare invariance 

In general, we are interested in d-dimcnsional backgrounds which have a four-dimensional Poincare 
symmetry (the restriction to four dimensions is unnecessary, but is the case most relevant for 
phenomenology) : 

ds^ = g^^dx" dx" = e'^^'^ab dx" dx^ + gij{z)dz' dz^ . (2.1) 
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Here, — (x", 2;*), where x", for a = 0, . . . , 3, are the usual coordinates of four-dimensional 
Minkowski space and — for i = 1, . . . , n, are the coordinates on the n — (d — 4)-dimensional 
transverse space. ^ We will assume that ^4(0) 0, so that the four-dimensional metric at the origin 
in the transverse space is canonically normalized. In this present work we will concentrate, for the 
most part, on a more restricted set of metrics which are conformally flat; that is of the form 

ds"^ = e-^(^) {rjab dx" dx^ - dz' dz') , (2.2) 

with a suitable choice of coordinates. Notice that when d = 5, an arbitrary metric of the form 
(2.1) is conformally flat, so in that case (2.2) is perfectly general. 

We will be interested, amongst other things, in smooth versions of the RS metric in five 
dimensions discussed in [1]. The metric is usually written in the form 

ds'' = e-^^'\ab dx" dx^ - dr^ (2.3) 

with A{r) = 2k\r\, but can be written in conformally-flat form with A{z) = A{r{z)) = 21og(A:|2;| + 
1). One simple way to introduce thick brane is by smoothing out RS ansatz. For example, one can 
make a substitution |r| — > /^"^ logcosh(/xr), although for practical purposes it is more convenient 
to smooth in the z-basis: 

\z\ — > — logcosh(/x^) . (2.4) 
A* 

Here /i is an independent parameter which determines the thickness of the brane. In the RS limit, 
II ^ k, we expect that any additional fields will be localized near z — 0, which will be assumed 
in much of the following discussion (we will comment on subtleties associated with smearing of 
the matter fields over the thickness of the brane later). On the other hand, for simplicity it is 
convenient to study the behaviour of the gravitational modes in a different limit fi ^ k, where 
A{z) depends on only one scale k, and we will do so in most of our examples. In this case the 
equivalent smoothing in conformally-fiat coordinates can be written as A{r) = 2 logcosh(A;r), or 
equivalently A{z) = log{k'^z^ + l)i [30]. This metric approaches the AdS form asymptotically for 
\z\ ^ 1/k. However, we will also consider metrics which are more general and not necessarily 
asymptotic to an AdS space. 

We first point out that smoothing of the RS solution can be performed without the addition of 

matter fields. Consider the five- dimensional case where the domain wall is generated by an explicit 
position dependent term in the gravity action, much in the spirit of the original Randall-Sundrum 
scenario. (We will later study gravity in the background of branes created by scalar fields, and 
find that the supersymmetric potential introduced in [27, 28] appears naturally in the solution of 
the field equations.) In the absence of fields other than gravity, we study the action 

S^-jd'x[^ {k-'R + A(r)) + ^)V{r)] . (2.5) 

^In our conventions the metric g^^, has signature (+,—,—,—,...). 



3 



Here g — \deig^y\^ ^'(4) = Idetgiabl is the determinant of the induced metric on the domain wall 
and — M~^, where is the fundamental Plank scale in five dimensions. The function V{r) 
will approximate a delta function which generates the domain wall, and A(r) is roughly constant 
away from the domain wall and corresponds to the bulk cosmological constant. This action gives 
rise to a stress tensor 

= \ (A(r) g,, + V{r) ga, 5« S^) . (2.6) 

We note that any stress tensor which is four dimensional Lorentz covariant can be decomposed in 
this way, and therefore is derivable from an action of the form (2.5). We are interested in actions 
which admit solutions of the form (2.3). The "bulk cosmological constant" and "brane tension" 
are then determined to be 

A(r) = -3K-M'(r)2 , V{r) = ?>K-^A"{r) . (2.7) 

For A{r) = 2k\r\, this reproduces the original Randall-Sundrum system [1]. By choice of a static 
background metric (2.3) the four-dimensional effective cosmological constant must vanish, and 
the fine tuning of cosmological constant and brane tension is replaced by (2.7). Indeed, one can 
integrate out the extra dimension and check that the action (2.5) vanishes for the background 
metric (2.3), which is equivalent to vanishing of the four-dimensional cosmological constant. 

We should comment that because this analysis does not depend on what type of field creates 
the domain walls, we can study solutions in which the domain walls have negative tension. For 
example, we can study smooth versions of the RS solution to the hierarchy problem [2], in which 
the transverse direction is compactified on an orbifold S^/Z2- Branes are stuck at the two fixed 
planes of the orbifold action, one of which has positive tension and the other negative tension. 
In the RS scenario the metric has the form (2.3) with A{r) = 2k\r\, —Vc < r < rc, where it is 
understood that r is periodic over an interval 2rc. In order to make the periodicity explicit, we 
can study a multiple cover of the circle, with 

A{r) = ^k[rc + (-1)" (2r - (2n + 1) Vc)] e{r - nrc) 9{{n + l)rc - r) , (2.8) 

neZ 

where 9{r) is the usual step function. The brane tension is related to A"{r), which has delta 
function singularities with positive coefficient for r = mod 2rc and negative coefficient for r — Tc 
mod 2rc. As before the solution can be smoothed by an appropriate smearing of the ^-functions; 
for example, 

e{x) e{x) = |(tanh(//x) + 1) , (2.9) 

where is, as described previously, a parameter which characterizes the thickness of the brane. 
An alternative smoothing procedure which may be useful for numerical calculation is a truncation 
of the Fourier expansion of the sawtooth: 

^ W ~ t - § E (2n\l)^ cos(7r(2n + l)r/r,) . (2.10) 
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We will also consider backgrounds that can be interpreted as three-branes embedded in spaces 
of dimension d > 5 and also backgrounds which can be interpreted as intersections of higher 
dimensional branes with a four-dimensional intersection. In the latter cases, the question for 
these kinds of backgrounds is whether four-dimensional gravity is localized on the intersection. 
For example, we will consider smooth versions of the multi-dimensional patched AdS space with 
metric [5] 

= , , Y^n I ^|X2 (^«fe dx^ - dz'dz') . (2.11) 

This metric represents n = c? — 4 intersecting [d — l)-branes, which mutually intersect in four 
dimensions. 



3. Universal Aspects 



3.1 Gravitational Fluctuations 



In this section we derive a universal equation for the effective four-dimensional gravitational fluc- 
tuations in the conformally-fiat backgrounds described in the last section. The more demanding 
case of the general background (2.1) will be discussed separately in Sec. 3.4. This means that we 
study fluctuations of the metric (2.2) of the form 

^ g-A(z) 1^^^^ ^ ^^^(^^ ^^a ^^b _ ^^i _ ^3_-^^ 

It will be convenient to define hf^^, to be a fiuctuation whose only non-zero components are hab- 
We will use the transverse traceless gauge for these fiuctuations, i.e. h^^ = and /?.^ = 0. We 
should note that there may be additional fluctuations not of the form (3.1) in transverse traceless 
gauge, but we will not comment on such modes here. 

Since the metric (2.2) is manifestly conformally flat, it is convenient to use the the general 
form of the Einstein tensor for metrics of the form g^^,^ = e.~^gni, (see for example [53]): 



G 1X1/ G fXV ~l~ 



d-2 



\V^A V.A + V. A - ~g^, ( V, - ^ VpA VM) , (3.2) 



where indices are raised and lowered with g^j_i, in this context. Using the form of the Einstein tensor 
for linear perturbations about flat spacetime [53], 

5G^, = d''di^,h,)p - \d^d, V - \d,d.K - h/^- (^'^" V - d'd, K) , (3.3) 
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the linearized perturbation of the Einstein tensor (3.2) becomes, 

+ ya^QM + 77^, {-dp h''^ d^A + \dp hi d^A - hP^ dpd^A) 



(3.4) 



Only the terms which are underlined are actually non-zero. The other terms either vanish due to 
the gauge conditions or they vanish because only has non-zero components hab and, moreover, 
is a function of the {z^} only. 

The next question concerns the variation of the stress tensor. For any action of the form (2.5) 



'-'-'-HI/ -'■ju "'Kl/ 1 



(3.5) 



which is automatically symmetric. Later we will show that this transformation property remains 
valid when the background is generated by scalar fields. From (3.5), and the unperturbed Einstein 
equation, we derive 



d^A + df,d^A) h^, + ^K-''{-dpdPA + ^ dpAd^h 



<i-3 



(3.6) 



Again only the underlined terms survive due to either the gauge conditions or the particular choice 
of background. Both the surviving terms of 5T^y cancel with two terms of 50^,1, in the graviton 
equation of motion 5G^y — k'^STh^ = 0, leaving 



kd''dph^, + ^d''Ad,h,,^0 . 



(3.7) 



If we redefine the metric perturbation so that its kinetic term has the canonical normalization, 
i.e. hp,jj — e*^^"^)"^/^ /i^,^, then the term linear in derivatives is removed: 



^^^d^AdpA - ^dPdpA 



(3.8) 



We should emphasize that indices are raised and lowered here using the fiat metric rj^i,. In addition, 
only the hab components of the fiuctuation h,j_i, are non- vanishing. 

Now we use the fact that di^dp = — V^, where = —r]"'''dadb and = df, and look 
for solutions of the form hab{x,z) = hab{x)ip{z) with \3jiab{x) = m'^hab{x), where m is the four- 
dimensional Kaluza-Klein mass of the fiuctuation. Then since A — A{z) only, we have 



-v'Mz) + 



16 



tp{z) = m^'ip{z) . 



(3.9) 



This has the form of a Schrodinger equation for the "wavefunction" ^l){z)., "energy" m and potential 



V{z) 



i.d-2f 
16 



V^A{z)-V,A{z)-^VlA{z) 



(3.10) 
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The fact that when expressed in terms of the variable hfj,i,, the graviton equations-of-motion 
(3.8) have no single derivative terms is equivalent to the fact that the action for these fluctuations 
has the form of a canonical kinetic term: 

S J d^xOph^^dPh^" , (3.11) 

where indices are contracted with the flat metric rifj_„. This can be seen by expanding the scalar 
curvature about the background (2.2) and keeping track of powers of the conformal factor in the 
metric. (If we had chosen coordinates in which the metric is not explicitly in the conformally-flat 
form, then the kinetic terms in the x and z directions would have had different conformal factors 
multiplying them.) The redefinition of the metric hn^, — e^'^~^^'^^^^/^/i^,^ then absorbs the conformal 
factor multiplying the kinetic terms and puts the action in the canonical form (3.11). In particular, 
for solutions in the form hab{x)ip{z), (3.11) includes the term 

J (Tzi^izf ■ J d^xdchab{x)d''h"^{x) , (3.12) 

from which we deduce that the appropriate inner- product for the "wavefunctions" ip{z) is the 
conventional quantum mechanical one. (Notice that this differs from the inner-product employed 
in [27].) 

In order to calculate the strength of the four-dimensional gravitational coupling, it is convenient 
to decompose the action into a four dimensional part and higher dimensional parts, before rescaling 
the graviton. As in [1,5], including the fundamental (d-dimensional) Planck scale M*, which is 
related to the coupling k via 

the action takes the form, 

S ~ Mt' f (Tz e-('^-2)^/2 . / d^x + • • • , (3.14) 



where ^(4) is the determinant of the four-dimensional metric for matter perturbations about flat 
spacetime, and R^'^^ is the four-dimensional curvature scalar created by those matter perturbations. 
This allows us to identify the four-dimensional Planck scale M4 via. 

Ml = j rz e-('^-2)^W/2 , (3.15) 



and determines the four dimensional gravitational coupling Gn M, 



-2 3 



^Note that the above relation rehes on our choice A{0) = 0. 
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3.2 The four-dimensional graviton 



The question of whether there is locahzed (four-dimensional) gravity supported in the vicinity 
of the brane now becomes contingent on properties of the quantum mechanical system described 
by the Schrodinger equation (3.9). In particular, in order to have an effective four-dimensional 
theory of gravity wc require that (3.9) admits a normalizable zero-energy ground state. To find 
this zero-energy state, we notice, generalizing the observation of [27,34] to higher dimensions, that 
the Schrodinger equation (3.9) can be rewritten as a supersymmetric quantum mechanics problem 
of the form, 

■QjP{z)^m'^{z) , (3.16) 

where the n- vector of supersymmetry charge is 

Q^V, + ^V,A, gt = _v, + i^V.yl . (3.17) 

Hence, the zero-energy wavefunction, annihilated by Q, is 

^j;,(z)^exp[-^A(z)] . (3.18) 

Notice that such a wavefunction always exists since (3.7) always admits the solution where hab = 
hab{x), only, and Hxhabix) — 0. Furthermore, since the "Hamiltonian" • Q is a positive definite 
Hermitian operator, there are no normalizable negative energy graviton modes, as required for 
stability of the gravitational background. 

The condition for having localized four-dimensional gravity is that ipo{z) is normalizable; in 
other words 

I (Tz exp [ - ^A{z)\ < oo . (3.19) 

Notice that normalizability of the ground state wavefunction is equivalent to the condition that 
the four-dimensional gravitational coupling (via (3.15)) be non- vanishing; indeed 

G.~"'^--;^.°'°r (3.20) 

(V'olV'o) 

This requires that A{z) ^ oo sufficiently fast as \z\ — > oo. Normalizability is intimately connected 
with the asymptotic behaviour of the potential of the Schrodinger equation (3.9). If V{z) > as 
1^1 — >• oo, then 'ipo{z) is always normalizable. On the contrary, if V{z) < as \z\ — >• oo, then ■4!o{z) 
is not normahzable and therefore is of no interest to us since it cannot describe locahzed four- 
dimensional gravity. The situation where V{z) = as \z\ —>■ oo, is, perhaps, the most interesting 
and we will focus on that case in most of the remainder of this paper. 

At this point, we make the obvious remark that in any scenario where the transverse space is 
asymptotically flat Euclidean space {A{z) — > constant, as \z\ oo) '^o(^) is non-normalizable and 
gravity cannot be localized. 



8 



3.3 Corrections to Newton's Law 



In order to have localized four-dimensional gravity, we also require that the other solutions of 
the Schrodinger equation (3.9), the KK modes, do not lead to unacceptably large corrections to 
Newton's Law in the four-dimensional theory. 

In any realistic brane scenario, the matter fields in the four-dimensional theory on the brane 
would be smeared over the width of the brane in the transverse space. Rather than deal with 
this complication, wc will, for simplicity, consider the gravitational potential between two point- 
like sources of mass Mi and M2 located at the origin, = 0, in the transverse space [1,5,34]. 
This assumption is justified in cases when the thickness of the brane is small compared with the 
bulk curvature. We expect that our conclusions will be at least qualitatively correct in a general 
case and present some supporting arguments in Sec. 6. In order to evaluate the correction to 
Newton's Law, we note that a discrete cigcnfunction (these are not present in the RS case) of 
the Schrodinger equation il'miz) of energy w? acts in four-dimensions like a field of mass m and 
consequently contributes a Yukawa-like correction to the four-dimensional gravitational potential 
between two masses Mi and M2: 

U(r)^G.^^Mr'^^^^^J<^? (3.21) 
r r 

where the wavefunction ipm{z) is normalized / (PziprnizY = 1.^ As long as m is large enough, 
this will be a small correction. The fact that ■0m(O) appears in (3.21) is due to our simplifying 
assumption that the sources are point-like and located at — 0. A more complete analysis would 
involve the effects of the overlap of the gravitational modes with the matter modes, and would 
correct the factors of iprn{0)- We will not have more to say about such corrections here. 

The correction from any continuum states ipm{z) is obtained by integrating over these states 
with the relevant density-of-states measure. For states which form a continuum in n-dimensions 
starting at mo the correction to Newton's Law is 

U(r) ~ + /" ^«^^„(0)= , (3.22) 

where the wavefunctions ipmiz) of the continuum are normalized as plane waves, i.e. to unity over 
a period at — > 00. The factor of m"~^ is just the n-dimensional plane wave continuum density 
of states (up to a constant angular factor) . Notice that in flat d-dimensional space there would be 
no normalizable zero-energy wavefunction xjjo{z) and the continuum would extend down to m = 
and be unsuppressed: ipm{0) = 1. In such a case U{r) ~ MiM2M^~'^r~"'~^ — MiM2M^~'^r^~'^, as 
expected for the gravitational potential in d- dimensions. 



^Notc that in our convention the zero-energy state 'ipoiz) is not unit normahzed; however, since we have chosen 
A{0) = we have ^o(O) = 1. 
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In the case of intersecting branes, there are continuum modes which are locahzed in a subset 
of the transverse dimensions. In other words, the wavef unctions behave as plane waves only in 
p < n of the transverse dimensions. In this case the corrections from this part of the continuum 
spectrum are of the form (3.22) with n replaced by p. 

For the case when V{z) > as \z\ oo, the excited states are clearly separated by a gap from 
the ground-state. Hence corrections to Newton's Law are exponentially suppressed as in (3.21). 
As we have already mentioned, the most interesting case is where the potential goes to zero at 
infinity. In this case there is a continuum of scattering states ipmi^) with eigenvalues > and 
so the behaviour of the soft modes at 2;* = is crucial for determining whether the corrections are 
small. 



3.4 Extension to non-conformally-flat backgrounds 

In this section, we briefly indicate how the preceding analysis of gravitational fluctuations extends 
to cases where the background is of the general non-conformally flat form (2.1). In order to derive 
an equation for metric fluctuations 

ds'^ = e-^^'\riab + habix, z)) dx" dx^ - gij{z)dz' dz^ , (3.23) 

we follow essentially the same steps as for the conformally-flat case detailed in Sec. 3.1. First of 
all, it is convenient to deflne the metric g^uiz) — e.^^^^ g ij,v{z) and apply the relation (3.2) in order 
to flnd the variation of the Einstein tensor G^i,. Rather than writing down all of the terms, as we 
did in Sec. 3.1, we will make immediate use of the following four facts: (i) d^h^^ — 0; (ii) = 0; 
(iii) (7^jy only depends on z; and (iv) the only non-vanishing components of the variation h^i^ are 
hab- By brute force one can show that the variation of the Einstein tensor G^i, is 

5G^, = -iV^'VpV = -l^'^cV - Ir'^^d.i^g'^djh^,) , (3.24) 

where, until further notice, indices are raised and lowered with Using (3.2), we can then write 
down the variation of the original Einstein tensor: 

+ - r'^%{V^g'^d,A) + ^d,Ad'A)h,, . 
Assuming that the variation of the stress-tensor is given by (3.5), we have 

5T^, = ^K-'i- g-'^%{y^~g'^d,A) + ^K-'diAd'A)h,, . (3.26) 
Hence, Einstein's equation gives 

-la-^^eV - ir'^%iV^~9''djh^u) + ^d.Ad'h^, = , (3.27) 
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which can be rewritten in terms of the original metric as 

^<9p(V^^'"^a„V)=0, (3.28) 

where now indices are raised and lowered with g^^,. In other words, the general equation for 
the fluctuations is simply the covariant scalar wave-equation. This has been noted in the case 
of one transverse dimension in [34] which also discusses its significance within the AdS/CFT 
correspondence. For a conformally-flat background (3.28) reduces to (3.7). 

For a fluctuation of the form hab{x, z) — (p{z)hab{x), with □^^(jb(x) = m^habix), we have 

-^di{^g^^ dMz)) = m'g'\{z) . (3.29) 

This describes the 2;-dependence of a mode with four-dimensional Kaluza-Klein mass m. Notice 
that when m — (3.29) always admits the solution (p{z) — constant; this will lead to the analogue 
of the zero-energy solution tj^oi^) in the conformally-flat case. The fluctuation equation (3.28) is 
derivable from the action 

J d'^xy/gd.h^^d^h^'' = J rzg'^\z)^^^{zf- j d^'x d,Kb{x)d'h''\x) + ■ ■ ■ . (3.30) 

From the z-integral above, inserting (p{z) — constant, we deduce the generalized expression for the 
zero-energy "wavefunction" i)Q{z): 

Uz) ^ [g'\z)^)]"\ (3.31) 
which reduces to (3.18) in the conformally-flat case. The normalizability condition is consequently 

j d"z^™(^)v^<oo . (3.32) 

The "wavefunction" is consequently 

i^{z)^^{z)[g'\z)^)]"\ (3.33) 
which satisfies a generalization of the Schrodinger equation (3.9): 



4. Gravity Localized on Thick Three-Branes 



In this section we consider in some detail the case when the 3-brane is embedded in a space 
with dimension d > b. We will, as per Sec. (2.2), restrict ourselves mainly to a conformally-flat 
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background which is, in addition, radially symmetric in the space transverse to the brane. In other 
words we shall consider metrics of the form 

ds^ = e-^(^) {r]ab dx" dx^ - dg" - o'd^f) , (4.1) 

where q is the radial coordinate in the transverse directions and Q are the angular coordinates 
on S*^"^. The function A{q), as indicated, only depends on the radial variable. We will briefly 
consider the more general radially symmetric background which is not necessarily conformally flat 
at the end of Sec. 4.4. 



4.1 Localization and decoupling in d — h 



We begin by discussing the case in d = 5, corresponding to the RS scenario, when the transverse 
space is one-dimensional. In this case, the metric (4.1) is 

ds" = e-^^'^ {riab dx" dx^ - dz^) . (4.2) 

The Schrodinger equation (3.9) is: 



The zero-energy state (3.18) is [28,30] 

Mz)^e^p[-lA{z)] . (4.4) 

For this to be normalizable, exp[— |A(2;)] must fall off faster than 1/z. As we explained in Sec. 3.2, 
the question of whether 4'o{^) is normalizable is intimately connected with the asymptotic be- 
haviour of the potential V{z). If V{z) > as \z\ oo, then ipo{z) is always normalizable. On the 
contrary, if V{z) < as \z\ —>■ oo, then ipo{z) is not normalizable and therefore is no interest to us 
since it cannot describe localized four-dimensional gravity. 

The fact that there are no bound-states with negative energy follows from the factorization 
(3.16): 



-^ + jA'(z]\^(z]^m'^{z) , (4.5) 



which has the form of supersymmetric quantum mechanics Q^Qip{z) — m'^il^{z), with Q = d/dz + 
jA'{z). Hence the zero-energy state (4.4), which satisfies the supersymmetric condition Q'ipo{z) — 
0, is the ground state, the bound-state of lowest energy. 

The borderline case, V{z) ^ as \z\ oo, is the most interesting case, and includes smoothed 
versions of the AdS scenario of [1]. For example, when the wavefunction (3.18) has a power-law 
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Figure 1: The Schrodinger potential (3.10) for the case A{z) = log(A;^z^ + 1). 

fall-off ip{z) ~ |^|~", then normalizability requires a > |. In this case the potential V{z) falls off 
as a{a + l)/-2^. Conversely, if we assume that the potential falls off as \z\^'^'^ the wavefunction 
ipo{z) falls off as e"'^/'^' \ Consequently, if /3 < 1, then the wavefunction is not normalizable. It is 
interesting to note that the borderline case, where the potential falls off just slowly enough to give 
rise to a normalizable bound-state, i.e. V{z) ~ includes the AdS scenario. Figure 1 shows 

the potential V{z) for a particular smoothing of the AdS case with A{z) — log{k'^z'^ + 1), discussed 
in Sec. 4.2. 

We now have to consider the question of whether or not the other modes in our problem 
decouple. The relevant effects of these modes on Newton's Law were established in Sec. 3.3. For 
the case when V{z) > 0, for \z\ oo, the excited states are clearly separated by a gap from the 
ground-state. Hence corrections to Newton's Law will be exponentially suppressed, as in (3.21). 
The most interesting case is where the potential goes to zero at infinity. In this case there is a 
continuum of scattering states ip,n{z) with eigenvalues > 0. Since the bottom of the continuum 
is at m = it is clear that decoupling is a delicate issue. From (3.22), given that gravity is 
localized so that i^o{z) is normalizable, decoupling would require that dme'^^^lipmi^)]"^ has no 
singularity at the bottom limit of integration. 

Before we attempt some rigorous analysis, let us consider the problem with some rather crude 
apparatus. If the continuum modes are to decouple we want the probability for continuum modes 
to tunnel into the central region of the potential to be vanishingly small for m — > 0. In order to get 
a feeling for what might be required, it is instructive to consider the WKB approximation for this 
tunneling probability. Consider a continuum mode with energy incident upon the potential 
from the right. The transition probability is, in the WKB approximation,^ 



T{m) ~ exp 



zi{m) 



-2 / dz ^/V{z) - m2 

zo(m) 



(4.6) 



■''Of course, the WKB approximation is not valid in the central region of the potential for states of small energy. 
Later we shall present a more rigorous analysis. 
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where < zo{m) < zi{m.) are the two points in the rightmost barrier region where V{z) — m^. In 
order for the soft continuum states to decouple we would need T'(O) = 0. Since zi{Qi) — oo, this 
can be achieved if the integral (4.6) has a divergence for large z; in other words, if V{z) falls off 
at least as slowly as 1/ z^. This matches precisely the condition on the normalizability of the state 
'iPq{z). This suggests that there is a natural decoupling of the continuum states precisely when the 
ground state '4>o{z) is normalizable. 

So the hypothesis that we want to estabhsh is that the continuum modes decouple (that is lead 
to small corrections to Newton's Law) precisely when i)Q{z) is normalizable. To this end, consider 
a potential of the form 

V(.) . , (4.7) 

Z'^ 

for large \z\. We shall find that the crossover from locahzation to de- localization occurs for some 
critical value of a. We shall not assume any particular form for the potential V(^z^ except, for 
simplicity, that it only depends on a single dimcnsionful scale A;, so that, for instance, the central 
region extends over a scale 1/A;; in other words (4.7) is valid for \z\ ^ Xjk. 

As we have discussed above, what we need to calculate in order to investigate the decoupling 
of the continuum modes is the limiting behaviour of ■0m(O) at small Let us consider four 
different regions in z (we will always be considering modes with energies m <^ k): (1) z <^ 1/k, 
(2) 1/A; < 2; < 1/m, (3) 1/A; < 2; ~ 1/m and (4) 1/m < 2;. In regions (2), (3) and (4), we must 
solve the Schrodinger equation with the tail of the potential (4.7): 



The solution is given by a linear combination of Bessel functions 

V'm(^) = a^z^l'^Y^^xi2{yriz) + 6^2;^/^ Ja+i/2(m2;) . (4.9) 
In region (4), where mz ^ 1, the Bessel functions become plane waves: 



■0m(-2) = am\l — sm(mz - f a - f ) + bm\ — cos(mz - f a - f ) . (4.10) 



In region (2), where mz <^ 1 (but A;2; ^ 1), the Bessel functions can be expanded in mz giving 
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«-i/2 V 2 y 

+ 



^m(^) = 

TT \mz J 



r(a + 3/2) V 2 



(4.11) 



^Recall from Sec. 3.3 that the wavefunctions ipmiz) must be normahzcd as plane waves for \z\ — !■ 00. 

''The corrections in the second square bracket, coming from the expansion of ^0+1/2(^2:) for smaU mz, are a 
power series in {mz)'^. The corrections in the first square bracket, coming from the expansion of Ya+i/2{'mz) are 
more compHcatcd since they depend on whether a + 1/2 is an integer, or not. However, the two terms indicated 
are the dominant terms for small mz. 
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In other words we can match the wavefunction in regions (2) and (4) by using the asymptotic 
behavior of the Bessel functions, and the exact form of the Bessel function then determines the 
behavior in the intermediate region, (3). 

Now we show how to match regions (1) and (2). In these regions mz <^ 1 and it is meaningful 
to solve the Schodinger equation as a series in m^. The first two terms are 

^m{z) = Mz) + m^{z) + • • • , (4.12) 

where '^q{z) is the suitably normahzed zero energy solution (3.18). The first correction (t){z) satisfies 
the inhomogeneous equation 

-^,+V{z)\^{z)^Uz) . (4.13) 

Now we can match (4.12) with z » 1/k, with the series in region (2) (4.11). In this region 
■0o(-2) ~ z~°', which matches the first term in (4.11) if ~ m"+^/^. The second and third terms 
in (4.11) then match the next term in the expansion (4.12) as long as h„i ~ ^-^+3/2 8 Ugij^g this 
matching procedure, we have determined the m-dcpcndence of a.^ and 6^. Since m is small, the 
dominant term in region (4) comes from the second term in (4.10) where the coefficient of the 
cosine goes like fe^m"^/^ ~ m~°+^. Hence to normalize the wavefunction il>ra{z) as a plane wave, 
we must multiple it by an overall factor of m"~^. The value of the wavefunction at ;2 = is then 
extracted from (4.12). To leading order in m, 

/ 777 \ a— 1 

M{Q) - (^) , (4.14) 

where the factor of k is dictated by dimensional analysis. Notice that in the AdS scenario of [1] , 
a = |, in which case we find ■0m(O)^ ~ m/k, in agreement with the exactly solved delta-function 
potential of [1] . 

Equation (4.14) establishes the asymptotic behaviour of the continuum modes at z ^ 1/A:, 
where k is characteristic decay width of the potential V{z). If the matter is smeared over the 
distances 1/k one needs to obtain more information about the wave-functions of the excited modes 
to study the corrections to the Newton's Law. It is important to note, however, that we obtained 
(4.14) without specifying details of the potential near z — Thus, we can easily introduce a 
potential such that matter is localized near z — yet the arguments leading to (4.14) remain 
valid. 

From (4.14), we find that the integral over the continuum modes in (3.22) is only non-singular 
at the bottom limit of integration if a > i. In this case the corrections are 

Uir) ~ G^—^ + ^3^,._, = ^^^(1 + Jk^i) ' (4.15) 



*A caveat is that the thu'd term in (4.11), coming from Ja+i/2(™-2): could match terms higher in the expansion 
(4.12); however, this would require a non-generic behaviour of the potential. 
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Figure 2: The shape of the potential (4.19). 



where C and C are dimensionless numbers. We are assuming that the metric only depends on one 
dimensionful parameter k, so that Gn ~ kM~^. Notice that when the potential falls off as (4.7) 
there are power-law corrections to Newton's Law with a universal exponent a determined simply 
by the long-range fall-off of the potential. We emphasize that the AdS case of [1] corresponds to 
taking a = |, which agrees with the exact calculation of the correction presented in [1]. 

Now wc sec that it is precisely when 'iPq{z) is normalizablc, i.e. when a > |, that the continuum 
states give a correction to Newton's Law which is suppressed relative to the leading term. 



4.2 Examples 



At this point it is probably worthwhile to consider some examples. To begin with, consider the 
class of conformally-flat five dimensional backgrounds for which 

>1(^) = f log(A;V + l) , (4.16) 

for some constants a. The case when a = | is particularly interesting because in this case we can 
easily transform back to the r coordinate in which case yl(r) = 2 log cosh(/cr). So when a = | 
the space given by (4.16) is asymptotically AdS and therefore represents a smoothing of the AdS 
example of [1], recently discussed in [30]. 



For general a, the potential in the Schrodinger equation and the ground-state corresponding 
to (4.16) are, respectively. 

In other words, this class of examples has precisely the asymptotic form of the potential discussed 
in the last section. The shape of the potential for a = | appears in Figure 1. 
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The example above has the advantage of being simple; however, since it only depends on one 
parameter, k, the thickness of the brane is comparable to the bulk curvature. Without any further 
assumptions, one expects that the matter fields living on the brane will be smeared over the distance 
1/k and such smearing may affect the estimates for the corrections to the Newton's law. In order 
to localize the matter fields near ^ = we can introduce a parameter characterizing thickness of 
the brane (as observed by the matter fields) in addition to the overall scale k. For example, in the 
RS scenario, we could replace |r| in A{r) = 2k\r\ by either rtanh(/ir) or log cosh (/xr), which 
both depend on a thickness paramater /i; however, in these cases we cannot calculate z = z{r) 
and hence V{z) in closed form. It is consequently more convenient to start in the 2;-coordinate 
basis with A{z) — ^log(A;|2;| + 1) (where the RS scenario has a = |) and smooth \z\ with either 
ztanh{ij,z) or ijL~^\ogcosh.{iJ,z). In the latter case 

A{z) = ^log(A;/x-Mogcosh(/x;s) + 1) , (4.18) 

from which we can easily calculate the potential 

^ ak/i^ A;(a + l)(cosh(2/iz) — 1) — 2/i — 2A; logcosh(/x;2) (4 19) 

2 cosh.^ {iJ,z)(^k log cosh.{ij,z) + iJ^Y 

which has the asymptotic form (4.7) for \z\ ^ max(//~^, k~^). The generic shape of this potential 
is illustrated in Figure 2 where we have introduced the parameters Vi = — V^(0), V2 = max[V{z)] 
and Zi, defined by V{zi) = 0. The depth of the central well is easily found to be Vi = akfi. In the 
limit jj, ^ k the brane is very thin and, for a = |, the potential approaches that of Randall and 
Sundrum [1]. In this limit for general a, V2 ~ A;^, independent of /x, while zi is asymptotically 

^i^^^log , (4.20) 

2// [a + l)k 

so the central well becomes more like a delta function as /i ':$> k. In this limit k controls the 
long-range fall-off of the potential, via the asymptotic form 

^ , a(a + l)k'^ , , 1 / , X 

So in the limit fi ^ k we would expect that matter fields arc localized near 2; = 0, in which case 
we can neglect the overlap of those fields with the gravity modes as corrections to Newton's law 
and the formulae of Sec. 3.3 will be vahd. The other hmit /J, <^ k, where the brane is much thicker 
than k, is also interesting. In this limit, both Vi and V2 scale like k/i, while zi asymptotes to 



2a 



-{kix)-^/'' , n^k . (4.22) 



In this limit we expect that the smearing of the matter fields over the transverse direction will 
become important and that the analysis of the corrections to Newton's Law as described in Sec. 3.3 
will require some modification. 
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Figure 3: The volcano box potential. 



4.3 Resonant modes 

The final issue that we mention regarding decoupling, is the possible existence of resonances. It 
can happen that for some particular energies, incident plane waves can resonate with the potential 
V{z) and consequently have a large value of ■0m (0). This possibility was also noticed in [30]. 

A useful toy model in which the Kaluza-Klein modes can be calculated exactly is the volcano 
box potential (Figure 3). Because the potential is exactly zero beyond the barrier the can be no 
bound-state at zero energy, but the depth and width of the well can be easily arranged such that 
there is a single bound-state, with a vanishing small energy rri^ < 0, and a continuum for > 0. 
The solution for a symmetric continuum wavefunction (orbifold boundary conditions, as in [2]) is: 




(4.23) 



where 



ki — -s/ m? -\- Vi 





(4.24) 



and the coefficients are given by 




c = 2k2k3 [^2 COS kizi ((e^'^^^i + e^''''')k3 cos k3Z2 + (e^^^^^i - e'^''''')k2 sin k3Z2) 

+ ki sin kiz^ ((e'*^^"i - e''^''')k3 cos k3Z2 + {e"^''' + e'^^'^)k2 sin A;3^2)] , 
d = ^'''t^X^^^^ [k2 sin kizi ((e^*^^^! + e^^^^^)k2 cos A;3^2 - {e^^^^' - e^^^^^)^ sin k3Z2) 
+ k2 cos kizi ( - (e^^^^^i - e^^^^^)k2 cos A;3^2 + (e^*^'^^ + e^^^^^)k3 sin k3Z2)] . 



(4.25b) 



(4.25a) 



(4.25c) 
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Resonant modes occur when the coefficient a of the growing exponential in the region zi < \z\ < Z2, 
vanishes, i.e. 

ki 

cos kiZi — — sin kiZi — . (4-26) 

We assume that < V2 <^ Vi. In order to make contact with smooth versions of the RS model we 
set ViZi = k and V2 = k'^, up to numerical coefficients. The smoothings can in general introduce 
other dimensionful quantities besides k; for instance, the width 2zi, or equivalently the depth Vi, 
of the well part of the potential. If we in addition take (ViZi)'^ <^ Vi (or Xl^/V2 <^ 1) then we can 
expand (4.26) in powers of (Vi-2i)^/Vi and m^/Vi to obtain. 

If there is a solution, then it is ~ 1^2 — (YiZiY ~ A;^. The spacing between resonances would be 
of order Ti'^/zf, so if there is a resonance below the barrier height V2 there is only one (for narrow 
wells). The contribution of the narrow resonance to Newton's law would be of the form 

U{r) . (4.28) 

Hence, the contribution of the resonance is negligible for r ^ Z2VV2 — rn? /m^. Notice that by 
(4.27) if ViZi > there will not be a resonance at all. In the RS case the potential is [1], 

1 "SP 

If we naively read off the coefficients 3/c and ^/c^ of the delta function and barrier height, respec- 
tively, then we find that there are no resonant modes in the volcano box approximation. However, 
this matching of the RS case to the volcano box approximation is too glib and a more careful 
analysis is needed. Hence, we find that if the well is deep, there is at most a single resonance 
at the order of, but below, the barrier height. We expect that this will be the case for general 
smoothings of the RS scenario, but in the absence of explicit solutions for the wavefunctions it 
is difficult to make definite predictions. In any case, because the resonances give rise to Yukawa 
type contributions to the gravitational potential, they are irrelevant below a certain scale which 
for smooth versions of the RS case is expected to be of order the fundamental Planck scale M*. 

For example, in Figure 4 we plot the "transmission coefficient" T — l/(c^ + d?) (which is not 
restricted to T < 1) into the well versus the energy of the KK mode w?. In this example, in units 
where the fundamental Planck scale M.^, = 1, we take Vi — 10^, V2 = 10, 2:1 = 1/Vi and Z2 = 10. 
We find a sharp resonance near = V2 — (Vi^i)^ = 9, as expected, of width (Am^)/m^ ~ 10~^. 
Away from the resonance T is smaller than O{10~^), except near V2, where T increases to .012 as 
a result of the weakening effect of the potential barrier at higher energies. 
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Figure 4: Resonance of transmission T into well with Kaluza-Klein mass m. 



4.4 Localization and decoupling in d > 5 



In this section, we consider the generahzation to the cases when the transverse space is more than 
five dimensional. In this case, it is convenient to use polar coordinates for z — {g, ft) and write 
ip{z) — R{g)Yi{fl), where Yi{D,) is an n-dimensional spherical harmonic, with 

l{l + d-6). 



Prom (3.9), we find that the radial function R{g) then satisfies 



-R"{g) - ^R'{g)+ 



i^A'igf - '-^A"{g) 
_ Mp)A'(^) + ^-^]R{g) = m'R{g) . 



(4.30) 



(4.31) 



The zero-energy state which potentially describes a four-dimensional graviton is from (3.18) 

Mq) = exp [ - ^A{g)] . (4.32) 



We now want to argue that, just as in the case in d — 5, the wavef unction ^00 (-2^) becomes 
normalizable precisely when the KK modes decouple. As in d = 5, the most delicate case is 
when ijjoi^) has a power-law fall-off ~ Normalizability requires that a > In this case 

the potential in (4.31) falls off as ~ g^'^. It is clear that only the s-wavc modes (/ = 0) are 
non- vanishing at the origin. Furthermore, we do not expect the matter fields of interest to us to 
transform under the SO(n) global symmetry and hence they would only couple at tree level to 
the s-wave modes. Finally, even for fields transforming non-trivially under SO(n) the sum over I 
presumably leads to a convergent series.^ For this reason we shall only consider the I — modes 
here. We now follow exactly same steps as in Sec. 4.1. Asymptotically for large g, in regions (2), (3) 
and (4), the radial function satisfies the equation, for I — 0, 

-R"ig) - ^R'ig) + ^^^^^R{g) = m^i?(^) . (4.33) 

^We thank Martin Gremm for raising the issue of the I > modes. 
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As previously, we can solve this equation in terms of Bessel functions: for Z = 

R{g) = a^^2-'^/V„+3_d/2(m^) + bmg'^-'^^^ Ja+3-d/2{mg) . (4.34) 

We now match the solution to that in region (1) using the same procedure that we followed in 
Sec. 4.1. In this case, we find after normalizing the wavefunctions as plane planes at ^ — > oo, 

/ Tfi \ a—d+i 

^m{0) ~ [j) , (4.35) 

Plugging this into the correction to Newton's Law, we find that the integral over the s-wave KK 
modes has no singularity if a > matching the condition that ipo is normalizable, and the 
correction is of the form 

^ M^M2 , C MM ^ M,M2 f^ , C ^ 

U{r) ~ Gn—^ + ^d_2^2,_2,+8 ^2a-d+5 = + ^kr^a-d+A ) ' (4-36) 

where C and C are some dimensionless numbers and we have assumed that the metric depends 
upon a single dimensionful parameter k so that Cat ~ M^'^k^"'^. 

Although we have only presented a detailed analysis for the conformally-fiat backgrounds, it 
is a simple matter, using the formulae of Sec. 3.4, to generalize to the arbitrary radially symmetric 
background having the form 

ds^ = e-^^^^riab dx'' - e-^(^) {dg'' + g'^dQ^) , (4.37) 

where the coordinates {g, Q} are polar coordinates in n = {d — 4)-dimensions. In the general 
case there is no canonical choice for the {z*} coordinates. It turns out that in order to analyze 
the localization of gravity it is not judicious to choose the {z'} to be the Cartesian coordinates 
associated to {g, Q}. On the contrary, we will choose the {z^} to be the Cartesian coordinates 
associated to polar coordinates {g,ft}, involving a new radial coordinate g — g{g), for which the 
metric (4.37) has the form 

ds'' = e-^(^) {riab dx" dx^ - dg") - e-^^^^g^dn^ , (4.38) 

where the functions A{g) and B{g) are related to A{g) and B{g) by the coordinate transformation 
on the radial coordinate. In the new set of coordinates {z^} the wavefunction is (3.31) 

M^) = exp [ - lA{~g) - = exp [ - ^A{~g)] , (4.39) 

where we have defined 

A{q) = ^^Mq) + . (4.40) 
To be completely explicit, the normalizability condition is 

J d'^ztPoizf = Vol(,S'^-^) J i'-^dg exp [ - ^A{g)\ < oo . (4.41) 



21 



The equation for the fluctuations (3.34) can be simphfled by separating the variables: il!{z) — 
R{g)Yi{Q,), which leads to a generahzation of the radial equation (4.31): 



-R"{g) - '-fR'{g)+ 



(d-2) 



4e 



U'(^) + 



R(g) = ni^Rig) . 



(4.42) 



When B(g) = A(g) the metric (4.38) is conformally-flat and the previous equation for the radial 
fluctuations for the conformally-flat case (4.31) is recovered. Our equation (4.42) matches that 
derived in [18] for the case when the transverse space is two dimensional. Notice that for s- 
waves the equation for the radial function R{g) (4.42) is identical to (4.31) with the replacement 
A{g) — > A{g). Moreover the expression for the zero-energy wavefunction 4'o{^) (4.39) is identical to 
that in the conformally-flat case with the same replacement A{g) — A{g). Consequently, we can 
use the same analysis as in the conformally-flat case to argue that gravity is localized, i.e. ipoi^) 
is normalized, when the s-wave continuum modes are decoupled. 



4.5 Thick three-branes from a scalar field 



In this section, we investigate whether the three-brane scenario that we have discussed in previous 
sections can actually be generated by gravity coupled to a single real scalar fleld. 

The flrst issue that we must verify is that our ansatz (3.5) for the behaviour of the stress tensor 
under gravitational fluctuations is, in fact, valid for a scalar fleld. The fact that it is valid follows 
from the dependence of the action of the scalar fleld on the metric. In general, the total action of 
gravity plus the scalar takes the form 

S = J d''x,/^[- K-^R + ld^<f)d''<f) - V(0)] . (4.43) 

The stress-tensor for the scalar fleld is 

T^. - ld^<P dA - Ig^u [^d^<P dpct> g'^^ - V(0)] . (4.44) 

For metric fluctuations h^j^i, which only have a non-vanishing component hghi since = (t){z) only, 
we have 

h^'^dA = . (4.45) 

Prom this and the form of the stress tensor (4.44), the variation of the stress tensor under met- 
ric fluctuations (3.5) follows immediately. In the coupled system, the analysis of fluctuations is 
naturally more involved. However, the fluctuations of the scalar are completely decoupled from 
the transverse traceless gravitational fluctuations as described in Sec. 3.1 and consequently all our 
previous conclusions regarding the localization of gravity are still valid. 
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Next we consider whether the three-branes can be formed from a single real scalar field. The 
system of scalar fields coupled to gravity was studied in the same context in [27-29,33]. Later, in 
Sec. 5.2, we will will prove that a single real scalar field cannot produce a three- brane in d > 5 
dimensions, so for the rest of this section we will take d = 5 and take the metric in the form (2.3) 
and take = 0(r) only. This ansatz, as we will see, completely determines the form of the scalar 
potential V(0) and the solution 4>{r). 

The graviton equation of motion is, as usual, G^j, = n^T^y, with T^y as in (4.44). Plugging in 
the ansatz (2.3) and = 0('"), there are two independent components of Einstein's equation: 

fi;20'(r)2 + 2/t2V[0(r)] + QA'{rf - QA"{r) = , (4.46a) 
K^cp'irf - 2K^V[(l){r)] - 6A'{rf = . (4.46b) 

It immediately follows that, 

K^V[(t){r)] = -3A'{rf + |A"(r) , (4.47a) 

fi;20'(^)2 = 3^/'^^^ _ (4_47b) 

Note that a solution only exists for 0(r) if A"{r) > 0. The scalar field equation following from the 
action (4.43) is 

-^M^a'" dMr)) + ^ = , (4.48) 



or, with our ansatz. 



-(f>"{r) + 2 A\r) (f>'{r) + = . (4.49) 



One can easily show that the scalar field equation is solved automatically by a solution of Einstein's 
equation. To see this we note that Einstein's equation imphes V ^T^^^ — 0, which itself implies the 
scalar field equation due to the general covariance of the scalar field action. ^° Hence, given the 
scalar potential V(r), any solution to Einstein's equation will automatically be a solution to the 
scalar field equations. Alternatively, any metric of the form (2.3) is a solution to the gravity and 
scalar cquations-of-motion if the scalar field has the form (4.47b) and the scalar potential is given 
by (4.47a). 

Furthermore, if 0(r) is a strictly monotonic function of r then we can we can implicitly define 

W[(t){r)\ = K-^A!{r) , (4.50) 

hence it follows from (4.47b) that 
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We are grateful to Shanta de Alwis for pointing this out. 



23 



Figure 5: The Schrodinger potential V{zi,Z2) for the metric (2.11) with k = l and \z\ z tanh(102). 

Then we can write the scalar potential as, 

V|^l = |(^)^-3m^)^ (4.52) 

We therefore see that the supersymmetric form of the scalar potential and BPS equations intro- 
duced in [27, 28] appear naturally in this approach. The ansatz (2.3) for the metric determines 
that the scalar potential can be written in the supersymmetric form (4.52). 

5. Gravity Localized on Thick Intersecting Branes 

In this section, we consider the possibility that gravity can be localized on the four- dimensional in- 
tersection of higher dimensional branes. In other words, we search for conformally-flat backgrounds 
(2.2) with d > 5, that we can interpret as a four- dimensional intersection of d — 4 {d — l)-branes. 

For example, we have in mind smoothings of the multi-dimensional AdS metric in (2.11). For 
instance, we could take \z\ ztanh(102;). For this particular smoothing in d = 6 with k = 1, 
Figure 5 shows the potential V{zi, Z2) appearing in the equation for the gravitational fluctuations 
(3.9). 

5.1 A solvable example 

In general, the Schrodinger equation (3.9) in the case of intersecting branes is fully n-dimensional 
and consequently rather complicated; however, in the case when 

n 

A{z)=Y,A^\^') . (5.1) 

i=l 
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Figure 6: Close-up of the potential near the brane intersection. 

it is separable by writing ip{z) = ip^^^z^) x ■ ■ ■ x case the solution represents n 

intersecting {d — 2)-branes where the intersection is four-dimensional. 

For each i, ip^'^\z^) satisfies a one-dimensional Schrodinger equation identical to (4.3) with 
A{z) A^'^z'). Hence, we can immediately draw on the results of Sec. 4 to establish the conditions 
under which gravity is localized on the intersection. Basically we require that exp [ — 
falls off faster than 1/z^ as \z^\ oo, for each i, so that the zero-energy state 

^oW=rf^(^^)x---x4")(z"), (5.2) 

is localized in all the transverse directions and therefore represents the four-dimensional graviton. 
In general, the spectrum consists of 2" different sectors, depending on whether each ip^^\z'^) is the 
normalizable wavefunction iI)q\z^) or a continuum wavefunction ipm{z^)- So for example, when 
there are 2 transverse directions there will be 4 sectors in the spectrum spanned by 

(1) 4')(.^)4^)(.^) , (2) v^S^)(.^)^S(^') , 

(3) ^S(^^)4'^(^^) , (4) ^S(^^)^S(^') • 

The first state (1) corresponds to the four-dimensional graviton. The set of continuum states (4) 
contributes to Newton's Law as in (3.22), 

6Mr) ~ Mt' / dm.dm, V^S (O)^V^S(O)^ , (5.3) 

Jo 

The set of continuum states (2) or (3) are localized in one direction and contribute as an integral 
over the single eigenvalue m2 and mi, respectively, 

A/T A/f r°° A/f A/f a—iT^^r 

S.^;U{r)-Ml-'k dm,^^ m^f + Ml-'k dm,^^ ^^(Of • (5-4) 

Jo f Jo f 
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5.2 Intersecting branes from scalar fields 

In this section, we consider whether the conformally-flat intersecting branes background can be 
generated from a single real scalar field. We argue that a single scalar field cannot produce a 
general intersecting brane background, but this leaves open the possibility that such backgrounds 
can be generated from more than one scalar field; for instance a single complex scalar as in [8]. 
In addition, we will also see, as promised in 4.5, that a single real scalar field cannot produce a 
conformally-fiat three-brane embedded in d > 5. 

We assume the metric has the form (2.2), and there is a single scalar field which generates the 
energy-momentum tensor, 

T,. = '^d,(P{z)d,c^{z) - lg^^{ld^c^{z) d,c^{z) - V[<P{z)]) . (5.5) 

The Einstein tensor can be read off of (3.2), and satisfies, for components transverse to the 
intersection, 

Gu = -Coo + ^ [l{diA{z)f + d^A{z)] . (5.6) 

The diagonal components of the stress tensor satisfy a similarly simple relation, 

Tu^-Too+{di<l)iz)Y , (5.7) 

so the particular combination Gu + Goo = i^^iTu + Too) of components of the Einstein equation 
relates derivatives of the scalar field to derivatives of the metric, independent of the form of the 
scalar potential V{4>): 

\K\d,<P{z)Y = ^[\d^A{z)diA{z) + dlA{z)\ . (5.8) 
The off-diagonal components of the Einstein equation are, 

Gij - K^Tij = ^ [\diA{z)djA{z) + didjA{z)\ - \K'di(t){z)dj(t){z) = , (5.9) 
or, using (5.8), 

[d,d,e^^^^'^f = [d^e^'-^y^] [dy^'^/'] . (5.10) 

This constrains the type of metric which can be obtained by a gravitating scalar field. In particular, 
the solution of (5.10) is, 

n 

A{z)=A{Y,(^iz') , (5.11) 

i=l 

for some constant coefficients Oj. By a redefinition of coordinates, this can always be recast in the 
form A(z) = A{z^), which preserves a (d— l)-dimensional Lorentz symmetry and hence can create 
only a single (d — 2)-brane in d dimensions. Therefore, it seems to be the case that additional fields 
are required to create intersecting branes in the presence of gravity. In addition, we conclude that 
a single real scalar field cannot produce a p-brane in d > p -|- 2. 
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6. Conclusions 



Wc have studied general features of gravity in domain wall backgrounds. For general domain wall 
type metrics we have found the conditions for there to be localized gravity on the domain wall 
or on the intersection of domain walls. It turns out that it is possible in generalizations of the 
Randall-Sundrum scenario for the graviton zero mode to be non-normalizable. It is precisely in that 
case that the effective gravitational coupling on the domain wall vanishes, and that the Kaluza- 
Klein modes become relevant at long distances. We have seen several illuminating examples in 
which smooth domain wall backgrounds are created by scalar fields, and in which the domain walls 
are created by non-dynamical sources in the absence of fields other than the graviton. Negative 
tension branes, such as appear in the RS solution to the hierarchy problem, can be studied in this 
formalism. We have studied intersecting domain walls in the absence of fields other than gravity, 
and we argued that a single scalar field is not sufficient to produce intersecting domain walls. We 
commented on resonant modes in the continuum of Kaluza-Klein modes, and argued that they are 
unimportant in smooth versions of the RS scenario. 

One issue that is worth commenting on is the fate of the Lykken-Randall scenario for the 
solution of the hierarchy problem [11], in the context of thick branes. In this scenario we consider 
again the general five-dimensional backgrounds (4.2) but now associate the matter fields that 
describe our world with another branc — the "TcV branc" — located at a point Zq in the transverse 
space. In order to create the necessary hierachy of scales we need Zq such that 

^2^g-2A(2o) ^ XeV^ , (6.1) 

where Mpi is the Planck mass. In the context of the backgrounds discussed in Sec. 4, this requires 
that Zq is much larger than k (since k ~ Mpi). In other words, the TeV brane sits at a place where 
we may approximate the potential V{z) by its asymptotic form (4.7). Since the relevant dynamics 
now takes place a,t z — Zq we have to re-assess the effects of the KK modes. In particular, the 
effective Newton's constant on the brane is now given by 

G^^Mm^. (6.2) 

This modifies the effects of the KK modes; in the case when the continuum starts at m = 0, the 
corrections from the KK modes (3.22) are modified to [11] 

In the above, the zero-energy wavefunction iJjo{z) must be normalized to unity {'4jq\iPq) = 1. As 
Zq increases then one might worry that the ratio 'ipmi^o) /i^oi^o) becomes larger and the effect of 
the KK modes is less suppressed and there would be large corrections to Newton's Law. Actually, 
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as explained in [11], this is not the case. Since zq hes in the tail of the potential (4.7), we can 
use the exact form of the solution in terms of Bessel functions (4.9) to assess the magnitude of 
the corrections. For small enough m the first term in (4.9) dominates, and one finds using the 
expansion (4.12) (taking into account the correct normalizations on the wavefunctions) 

— - ~ ^ oc m , (6.4) 
V'o(^o) V'o(O) 

and so the suppression is unchanged from the situation where matter fields arc located at z = 0. 
Small enough m in this context means that we can approximate ipmi^o) by the first term in (4.12). 

The discussion above, where the matter fields are located ai z — Zq but are still point-like in the 

fifth dimension, also suggests that in a more realistic scenario, where the matter fields are smeared 
out in the fifth dimension, the corrections to Newton's law will also be similarly suppressed and 
our previous insistence that the matter sources were point-like and located at z = was not overly 
simplistic. 
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